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Abstract 

The consensus model of Krause and Hegselmann can be naturally 
extended to the case in which opinions are integer instead of real 
numbers. Our algorithm is much faster than the original version and 
thus more suitable for applications. For the case of a society in which 
everybody can talk to everybody else, we find that the chance to reach 
consensus is much higher as compared to other models; if the number 
of possible opinions Q < 7, in fact, consensus is always reached, which 
might explain the stability of political coalitions with more than three 
or four parties. For Q > 7 the number S of surviving opinions is 
approximately the same independently of the size A'^ of the population, 
as long as Q < A^. We considered as well the more realistic case of a 
society structured like a Barabasi- Albert network; here the consensus 
threshold depends on the outdegree of the nodes and we find a simple 
scaling law for S, as observed for the discretized Deffuant model. 
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1 Introduction 

Can statistical mechanics help to describe opinion dynamics? The last few 
years have witnessed several attempts in this direction [H El El 111 El Ej and 
Monte Carlo simulations have become an important part of sociophysics [Zj, 
enlarging the new field of interdisciplinary applications of statistical physics 
jH]. The starting point is represented by a random distribution of opinions, 
which can be integer or real numbers, among a group of persons, or agents. 
Next, some simple dynamical mechanism is introduced so that, due to in- 
teractions between the agents, the opinion of each agent changes over the 



1 



time until, at some stage, a configuration is attained where the opinions are 
no longer modified by the dynamics and then remain the same. It is then 
interesting to study the possible stable opinion configurations at the end 
of the process to see, for instance, if it is possible that all agents stick to 
one and the same opinion ("consensus") or whether a polarization around 
fewer dominant opinions takes place. The most studied consensus models 
are those of Sznajd jj], Galam [H], Deffuant et al. j2| and the one of Krause 
and Hegselmann (KH) [Hj. The crucial differences between the models are 
the dynamical mechanisms to update the opinions, but they differ as well in 
other important aspects. For instance, in the original versions of the models 
of Deffuant and KH, the opinion variable is a real number between and 1, 
whereas in the models of Sznajd and Galam it is an integer (this possibility 
was also studied for Deffuant in JU])- Besides, there are as well differences 
in the way the topology of the society is conceived. In the Sznajd model 
the agents sit on a lattice and can have opinion-affecting interactions only 
with their lattice neighbours; in the models of Deffuant and KH, instead, one 
assumes a society in which each agent has the same probability to interact 
with everybody, although recently scale free network topologies have also 
been considered for Deffuant |J1 US] . Among the above-mentioned consensus 
models, those of Sznajd and Deffuant are meanwhile quite well-known, as 
the convergence to the final configuration is relatively quick, which allows to 
simulate populations with millions of agents 0. 

In this paper we focus instead on the KH model, which has not been 
investigated by many people so far. In its original version [3] one introduces 
a real parameter e, called confidence bound. At every Monte Carlo step, the 
randomly selected agent i with opinion Si takes the average of the opinions 
of those agents j such that \si — Sj\ < e. This averaging process makes the 
algorithm very time-consuming compared, for instance, to Deffuant, and is 
essentially the reason why most people of the computational sociophysics 
community do not find it attractive. However, the real interest behind con- 
sensus models is whether they are able to describe real situations, and it 
is not said that the faster the algorithms the better they are; the Sznajd 
model could effectively simulate the distribution of votes among candidates 
in Brazilian and Indian elections ^] . 

We study here a modified version of the KH algorithm, where opinions 
take integer values so that each individual has a finite number of possible 
choices. This is most often the case in real life; thinking for instance about 
elections, the voters have a limited number of possible parties and/or candi- 
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dates among which to choose. We have used several values for the number 
of opinions Q and checked how many different opinions survive in the final 
configuration. Initially we have assumed a society in which every agent in- 
teracts with all the others. This is however quite an unrealistic situation; 
therefore we have as well checked what happens if the personal relationships 
within the society form a scale free network, with few people having lots of 
friends and many having just a few. To build the network we adopted the 
popular "rich get richer" strategy proposed by Barabasi and Albert 



2 The model 

Our opinions can take the values 1, 2, 3, ...Q. As far as the confidence bound 
e is concerned, we will assume in this paper that the agents are influenced 
only by the individuals whose opinions differ by at most one unit from theirs. 
This corresponds to the case e = 1/Q m the original KH model. Actually, 
in order to go smoothly to the continuous limit one should introduce an- 
other parameter L, which is the discrete confidence bound, i.e. the maximal 
distance between compatible opinions (for us L = 1), and take the limit 
L, Q — > oo by keeping e = L/Q fixed. 

The algorithm starts by randomly distributing the opinions among the 
agents. The use of integer-valued opinions spoils the original KH concept of 
" average of compatible opinions" , because such average in most cases would 
not be an integer. 

A possible way out is to re-interpret the spirit of the original KH model 
in a probabilistic fashion. Suppose we want to update the status of agent 
i, which has opinion k. The number of agents with compatible opinions are 
Uk-ii and n^+i (respectively for opinions k — 1, k and A; + 1). If the total 
number of compatible individuals is n = n^-i +nk + nk+i, we say that agent 
i takes opinion k — l,koTk + l with probability Pk-i = rik-i/n, pk = rik/n 
and pk+i = Uk+i/n, respectively. This is to our mind a natural extension 
of the KH model to discrete opinions, and is the version we have used here. 
The status of the agents is updated sequentially, in an ordered sweep over 
the whole population; the program stops if no agent changed opinion during 
an iteration. 

The fact that the opinions are discretized allowed us to speed up the 
algorithm compared to the continuous case. In the latter the time to complete 
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Figure 1: Number of evolution steps necessary to have convergence to the 
final configuration as a function of Q, for the standard KH model with con- 
tinuous opinions and for our discretized version. The number of agents is 
1000 and for each value of Q we averaged over 1000 realizations. 



an iteration goes as N"^ {N is the size of the population), because for each 
agent to update one needs to make a sweep over the whole population to 
look for compatible individuals and calculate the average of their opinions. 
Here only the probabilities pk matter, so we keep an array where the opinion 
histogram Uk is stored {k = 1,2,..,Q). When we update the agent i with 
opinion k, from the histogram we derive directly the probabilities for the 
agent to take its next opinion. Suppose that agent i takes opinion j, what 
we need to do is to increase rij and to decrease Uk by one unit, to get the new 
opinion distribution, after that we can proceed to update a new agent. In 
this way we avoid to count each time the number of compatible individuals, 
which is very time consuming, and the time needed to complete the iteration 
goes as (2L + 1)A^ (for each agent one needs to make a sweep over the 2L + 1 
compatible opinion channels to get the probabilities, in our case 2L-|- 1 = 3). 
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Figure 2: Fraction of final configurations in which no consensus is reached, 
as a function of the number of agents, for Q = 7, 8. 

With our algorithm systems with millions of agents, unreachable by standard 
KH, can be simulated (it took us less than six hours to simulate one million 
agents on a PC). 

We found that the convergence to a stable configuration is much slower 
than in the original KH model. In Fig. ^ we plot the average number 
of sweeps necessary for convergence in the two cases, as a function of Q 
(we use the relation e = l/Q to make a correspondence between the two 
models). We see that, except for very high values of Q, the number of 
evolution steps for our algorithm is an order of magnitude higher than for the 
continuous model. This is most likely due to the stochastic character of our 
procedure; the opinion distributions vary more slowly if we allow jumps from 
one opinion channel to the neighbouring ones with some probability instead 
of sistematically shifting every agent to the average channel. Moreover, the 
standard deviation of the average evolution time is much larger in the discrete 
than in the continuous model, which hints to the presence of wild fluctuations. 
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Figure 3: Number of surviving opinions as a function of Q for a society 
where everybody interacts with everybody. We averaged over 1000 reahza- 
tions. 

3 Results 



The main result of our simulations is the existence of a threshold Qc, such 
that, for Q < Qc, consensus is always reached. This is true for both social 
topologies we have considered. For a society where each agent has relation- 
ships with all others, we find that Qc = 7 as we can see in Fig. El Here 
we plot the probability of having polarization as a function of the number 
of agents, for Q = 7,8. By polarization we mean that more than just a 
single opinion survive in the final configuration. The probability is given by 
the fraction of configurations with polarization. For Q = 7 this probability 
decreases strongly with and for = 10000 all samples presented a single 
final opinion. For Q = 8, instead, the probability for polarization is basically 
one for = 10000. We remark that our threshold is higher than in the Sz- 
najd model, where Qc lies between 3 and 4 and in the discretized version 
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Figure 4: Threshold for complete consensus as a function of the outdegree 
m of the network. 

of Deffuant^, where Qc = 2. This shows that the dynamics of the KH model 
is the most suitable to explain how competing factions can find an agreement 
and to justify the stability of political coalitions with several parties like in 
Italy. 

In Fig. ini we show how the number of surviving opinions S varies with 
Q, for several A^'s. We see that, as long as Q<^N, so that no finite size 
effects take over, S is approximately the same independently of the number 
of agents; the result holds for the Deffuant model as well [T^ . 

Let us now check what happens if we put the agents on a scale free 
network a la Barabasi-Albert. To build the network we must specify the 
outdegree m of the nodes, i.e. the number of edges which originate from 
a node. The procedure is dynamic; one starts from m nodes which are all 
connected to each other and adds further N — m nodes one at a time. When 

^In the original continuous versions the transition value of the confidence bound for 
consensus is ec ^ 0.4 for Deffuant and ec ^ 0.21 (our estimate) for KH. 
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Figure 5: As Fig. El but for a society where personal relationships are the 
edges of a Barabasi-Albert network. We averaged over 1000 realizations. 



a new node is added, it selects m of the preexisting nodes as neighbours, so 
that the probability to get linked to a node is proportional to the number 
of its neighbours. In all networks created in this way the number of agents 
with degree k, i.e. having k neighbours, is proportional to 1/k^ for k large, 
independently of m. 

In our simulations we took the network as undirected, so that commu- 
nication between two neighbouring agents can take place in both directions. 
We find again that there is a threshold Qc for the system to evolve to com- 
plete consensus. Interestingly, Qc depends on the outdegree m, as shown in 
Fig. 13 In fact, the final stable configurations are those in which each agent 
is surrounded only by agents which share its opinion or are uncompatible, a 
solution of a special graph colouring problem; only in this case each agent 
will maintain its opinion in the future with probability one. If m is small the 
average degree is small and it is easier to reach such configurations even when 
just a few opinions are available. Looking at Fig. |3J we see that Qc = 3, 4 
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Figure 6: Scaled plot of the same data as in Fig. |31 Apart from deviations for 
very small Q, where 5* is the same independently of A^, a reasonable scaling 
is observed. The two straight lines represent the two situations where few 
opinions survive (horizonal) and every agent keeps its own opinion (skew). 

for m = 1,2, respectively^. No polarization is possible for Q = 2 because, 
the network being connected, there would be at least two clusters sharing a 
border. As each agent has at least m neighbours by construction, for large m 
we expected to reach the threshold Qc = 7 that we have found in the case in 
which everybody is connected to everybody; this is indeed true for m larger 
than about 40. 

The analysis of the number of surviving opinions S is instead relatively 
independent of the outdegree m. We chose m = 3, in order to make compar- 
isons with corresponding results for the discretized Deffuant model JU] • Fig- 
El shows S" as a function of Q for different population sizes. The pattern looks 
very similar to the one observed in [10] for Deffuant: for Q not too small, S 

^As a matter of fact, in the special case to = 1 we find that consensus is not complete 
for Q — 3, but is reached in about 80% of the cases. 
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equals Q and only when Q gets close to A^, finite size effects take over and 
S converges towards A^. In the latter case, the simulation stops very early 
because most agents have different opinions and therefore the chance for an 
agent to change its mind is small. In ^lOj a simple scaling behaviour of 5" 
with Q and A*" was observed. The ansatz was 

S={Q~ l)f{Q/N)- fix ^ 0) = 1, fix ^ oo) = 1/x (1) 

In Fig. ini we rescaled the data of Fig. El according to the ansatz of Eq. 
n The "teeth" below the horizontal line refer to small values of Q, and we 
know that here S equals one or is close to one, independently of N; otherwise 
the scahng is good. 



4 Conclusions 



We have studied an extension of the Krause-Hegselmann consensus model to 
integer-valued opinions, both when all agents talk to each other and when 
they sit on the nodes of an undirected Barabasi- Albert network. We assumed 
that only agents with opinions differing by at most one unit can influence 
each other. A non-trivial implementation is necessarily probabilistic and 
many more iterations are required for convergence compared to the standard 
model. On the other hand, our algorithm is much faster than the continuous 
version and therefore larger population sizes can be explored. In a society 
where each agent can interact with all others, when no more than seven 
different opinions/positions are possible, the system always evolves towards 
consensus if the size A^ of the population is large enough (of the order of 10^ 
agents or more). On the other hand, on a network-structured society, the 
threshold Qc depends on the minimal number of friends an agent can have; 
if this minimum is less than 30, consensus is more difficult and Qc can lower 
up to 3. This versatility of the model makes it more suitable than others in 
order to explain how consensus can be reached in a variety of situations. In 
Italy, for instance, the 80's were the years of the so-called "pentapartito", a 
government's coalition of five parties. The stability of this coalition could 
be justified neither by the Sznajd model nor by Deffuant, but it is natural 
in our model (although the concept of political stability in a country where 
fifty-seven governments alternated in fifty-eight years is questionable!). 



10 



In a fully connected society, the number of final opinions S is an intensive 
quantity, i.e. independent of N for large N. 

On a network S grows with N if Q and increase so that the ratio Q/N 
is constant, but is intensive if Q is kept fixed when N —>■ oo. When Q is 
not too small 5* is well described by the same simple scaling function that 
reproduces the data for the discretized Deffuant model [TD|. Moreover, we 
find that S = Q for Q of the order of ten or larger, so in a realistic society 
with more than ten different opinions/positions, the dynamics of the KH 
model is unable to suppress any of them. 
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